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Abstract 



When a sizable, hard-sphere-like impurity is placed in a trapped Bose- 
Einstein condensate, it will create a hole in the condensate wave function. 
Since it costs more energy to drill a hole at the wave function's maximum 
(near the trap center) than at its minimum (near the trap edge), the impurity 
will be expelled from the condensate. This suggests that Ceo molecules freely 
falling through a condensate of atomic hydrogen may be deflected by a radial 
acceleration comparable to g. The scattered fullerenes would then form a 
characteristic banded density pattern outside of the trap. The detection of 
these deflected molecules can be used to identify the onset of Bose-Einstein 
condensation. 
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The Bose-Einstein condensation (BEC) of particles is the macroscopic occupation of a 
single quantum state. When the effect of interaction is weak, this occupation magnifies 
the ground state energy by a factor iV equal to the number of condensed particles. For 
macroscopic N m 10 24 , this extraordinarily large factor can be exploited to magnify minute 
changes in the ground state. Thus to the extent that the environment affects the condensate's 
ground state, BEC is potentially an exceedingly sensitive probe of its environment. 

In this work, we consider the effect of placing a hard sphere impurity into the condensate, 
which forces a hole in the condensate wave function. For translationally invariant systems, 
such as bulk liquid helium, it is immaterial where this hole is located. However, for atomic 
trap experiments |3]-|1J, because the condensate wave function is inhomogeneous, there is 
a position-dependent energy penalty. This energy penalty is ordinarily too small to be 
detected. We show here that such a minute energy difference can be magnified by BEC into 
a macroscopic observable phenomenon. 

This effect is generic, depending only on the inhomogeneity of the ground state wave 
function and its macroscopic occupation. It is sufficient to examine the problem in its 
simplest conceptual context. Consider a non-interacting Bose gas of mass m confined to a 
spherical cavity (an infinite square well) of radius b. An impurity, regarded as a hard sphere 
of radius a, is placed at a distance d from the cavity center. The situation is illustrated in 
Fig.l. The cavity traps the Bose gas, but it may or may not also confine the impurity. The 
impurity can, but need not be microscopic. The effect is actually most pronounced when 
a is comparable to b. Classically, there is no interaction between the hard sphere and the 
cavity. However, when the cavity is filled with a Bose gas, the ground state energy is altered 
by the presence of the impurity and there is a quantum mechanically induced interaction 
between the two. 

We will assume that the impurity is sufficiently massive that the Born-Oppenheimer 
approximation is adequate. Unless noted, distances are measured in units of b and energies 
in units of The effective impurity potential is then V(d) = NE , where E is the 

ground state energy of a single particle confined in a unit spherical cavity with an off-center 
hole: 

- V 2 ^(r)= EiI)(t), 
i/j(r) = at |r| = 1 and |r — d| = a. (1) 

When d = 0, the ground state energy and wave function are known analytically: E = 
7r 2 / (1— a) 2 , ipo(r) = - sin[7r(r— a)/(l— a)]. When d > 0, spherical symmetry is broken and the 
problem is non-trivial. While this can still be solved by some brute-force numerical method, 
the recent suggested use of a conformal transformation |J is more elegant. Since cylindrical 
symmetry remains intact, the wave function can be written as ip(r) = p^ 1//2 e* m£ 'i?(p, z), with 



dz 2 dp 2 p 2 
Let Z = z + ip, Z = z + ip, the transformation [| 

s. Z-c 



Q2 f)2 m 2 _ 1. 

' ' R(p,z)=ER(p,z). (2) 
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with c real, will map the unit circle z 2 +p 2 = 1 to the unit circle z 2 + p 2 = 1 and the off-center 
circle [z — d) 2 + p 2 = a 2 to the on-center circle z 2 + p 2 = a 2 with a = (d + a — c)/(l — cd — ca), 



provided that c = (1 + d 2 — a 2 — y (1 + d 2 — a 2 ) 2 — Ad 2 ) / (2d). Under this transformation, 
both the two dimensional Laplacian and the 1/p 2 term in ([J) are transformed with the same 
multiplicative factor J(r) = [ (1 + zc) 2 + (pc) 2 ] 2 /(l — c 2 ) 2 , 



d d m 



i ™2 _ I. - d 2 d 2 m 2 - - 



+ 



J(r) 



-+ 4 



dz 2 dp 2 p 2 J L $2 2 <9p 2 p 2 

Thus in terms of the transformed coordinates, t/>(?) satisfies 

J(f)(-V 2 )^(f) = EV(f), 
■0(f) = at |f| = 1 and |f | = a. (4) 

This equation can be solved by expanding |)],[7j 

^(f) = E^rV,(f), (5) 

and diagonalizing the resulting matrix 

My = kiifalJl^kj, (6) 

where 4>i(r) and kf are the complete set of eigenfunctions and eigenvalues of (§]) with J = 1. 

The resulting ground state energies from this transformed Hamiltonian are shown as 
black dots in Fig. 2. eo = vr 2 is the ground state energy of the cavity without the impurity. 
The matrix (H) is diagonalized using 100 basis states. For all impurity sizes, the energy 
is lowered as the impurity moves off-center. The intuitive reason for this is stated in the 
abstract. When the impurity approaches the cavity surface, d — > 1 — a, a — *■ 1, the basis 
states (pi are tightly squeezed and the method converges slowly. This poor convergence shows 
up in Fig. 2 as an unphysical up-turn in energy near d ~ 1 — a for a = 1/3 and a = 1/5. 
At the point of contact, the basis states are squeezed out of existence and the method fails. 
Thus if the impurity is not confined by the cavity, this method cannot be used for d > 1 — a. 
For the latter situation, the following variational trial function gives an excellent account of 
the off-center energy: 

^(r) = (l-^j)(l-r)e-^ a (7) 

The variational parameter a jumps quickly from zero to some a-dependent value for cPO.l. 
The resulting energies for a = 1/3 and a = 1/5, as computed by the Monte Carlo method, 
are shown as dash lines in Fig.2. The energy is continuous and smooth across the contact 
point and is only slightly above the diagonalization results. 

When the impurity is microscopic, a -C 0.1, both of the above methods are inefficient 
and the energy shift cannot be accurately computed. In this case, we follow the original 
suggestion of Huang and Yang and replace the impurity boundary condition by a pseu- 
dopotential. Thus instead of (]l|), we solve, 
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(-V 2 + Aixa5 3 (r - d ) d \f_ d \ \ r ~ d|)^(r) = Eip(r), 

ip{r) = at |r| = 1, (8) 



perturbatively in powers of a. The result, 



E, 



— = 1 + 2aj 2 (rtd) + (2a) 2 j 3 (%d) C os(nd) - -Mitd) , (9) 
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is shown as indicated in Fig. 2. The first order term is proportional to the cavity ground state 
density, which makes precise the notion that the energy penalty is greatest at the maximum 
of the ground state wave function. Note however that, for a > 1/5, the actual energy is much 
steeper than the first or second order results. The second order term above is from summing 
over only the s-wave intermediate states. As pointed out by Huang and Yang ||, such a 
second order energy shift can be positive because the pseudopotential is not hermitian and 
a singular negative term is removed by the derivative operation on the wave function. Since 
we have summed only over the s-wave intermediate states, the negative singular term in i/j^ , 
which should be proportional to l/\r — d|, started out as proportional to 1/r. Accordingly, 
we removed this contribution to the second order energy. This s-wave-only second order 
energy, while in principle incomplete, is exactly correct in the limit of d —>■ 0. As shown 
in Fig. 2, its inclusion yielded excellent agreement with numerical results at all values of d 
for a < 1/10. In comparison, the first order result is good only for a < 1/20. The effect of 
interaction among the condensed atoms can be incorporated via the Gross-Pitaevskii mean- 
field. The result is a slight energy shift || in accordance with the sign of the scattering 
length. For a condensate of atomic hydrogen, this effect is particularly negligible because of 
hydrogen's anomalously small s-wave scattering length [|J. 

An impurity of radius a and mass M will therefore be expelled from the condensate with 
radial acceleration 

ar ~ M2mb^ b dd r 1 > 

Since the fall-off in energy is increasingly steep as a approaches 1, the expulsion force is 
greater for larger impurities. If the impurity is microscopic, with a C 0.1, the perturbative 
result (D can be used. For an order-of-magnitude estimate we can keep only the first order 
term and replace the derivative by its average value, 

We have restored all factors of fi/m and b so that the overall dependence on b is clear. Note 
that this radial acceleration is proportional to the condensate density ~ N/b 3 . Currently, 
all experimental traps are highly cylindrical. For these cases, we can make an equal volume 
replacement: b 3 — > jb 2 L, where L is the length of the cylinder, 



An 2 h 2 N 



ar ^3Mm^LU> (12) 
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This corresponds to an impurity expelled radially from the middle of the cylinder. The 
numerical value |7r 2 = 13.16 is in good agreement with the exact cylindrical value of 
4/ Jl(zi) = 14.84, where z\ is the first zero of Jq{z). 

If the impurity is microscopic, the expulsion will be most pronounced for an impurity 
with a large radius-to-mass ratio, i.e., a hollow molecule. For Ceo, we have M ~ 10 -24 kg 
and asi4x 10~ 10 m. The latest reported Bose condensate in atomic hydrogen by Fried et al 
0] is ideal with optimal h/m = 6.65 x 10 -42 J-m 2 and the largest reported N « 10 9 . This 
condensate is highly elongated with cylindrical radius b = 7.5/im and length L = 5 mm. If 
C 60 molecules were sprinkled over this condensate, they would be deflected by an average 
radial acceleration of 

a r « 17m/s 2 « 2g. (13) 

For Na and Rb condensates, the effect would be reduced by a factor of 1/23 and 1/87 
respectively, modulo differences in condensate density. Since the Li condensate will collapse 
beyond a critical N c 10 3 , it will not be a favored condensate to observe this effect. 
To model the fullerene's additional surface attraction for hydrogen one should properly 
interpret a as the exact scattering length. Thus even if the surface attraction were to reduce 
the equivalent hard sphere radius by a factor of 2, the radial acceleration would still be 
comparable to g. 

The actual radial acceleration can be computed from the cylindrical energy shift 

h 2 4a J*( Zl d/b) 
E °- eo= m¥L-jK^T- (U) 

Fig. 3 shows the trajectories of C 60 molecules released uniformly at a height of 6b = 45fim 
above the cylindrical center and freely falling through the hydrogen condensate. They are 
narrowly focused near the trap's edge at ~ 15// m below the trap center. Below this point the 
beam diverges but is still bounded by two caustics. The density distribution at 66 = 45/im 
below the cylindrical center is shown in Fig.4. Most of the fullerenes are expelled out of 
the trap forming a sharply defined band just outside of the trap. The surprise is that the 
distribution is not diffuse, but is bounded by two high density caustics. This characteristic 
distribution is also generic. A similar band is also found using harmonic oscillator wave 
functions. For weaker repulsions (such as that exerted by the Na condensate), the focal 
plane and the caustics would form further down from the trap, but the banded density 
distribution would remain unchanged. 

We emphasize that this deflection is a purely quantum mechanical phenomenon. Clas- 
sically, if Cgo molecules were raining down on such a low density hydrogen vapor, only one 
out of 100 would encounter a single hydrogen atom. 

Since this repulsion is effective only when the atoms are condensed in the ground state, 
the detection of these deflected fullerenes can be used to identify the onset of Bose-Einstein 
condensation. Finally, if the impurity has a negative scattering length with respect to the 
condensed atoms, the impurity will be attracted and centrally focused by the condensate. In 
either case, a uniform distribution of impurities will become patterned after passing through 
a condensate. 
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FIGURES 





FIG. 2. The ground state energy of a spherical cavity with an off-center hard-sphere, a and d 
are measured in units of b. eo is the ground state energy without the hard-sphere. 
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FIG. 3. Classical trajectories of Ceo molecules falling through a cylindrical condensate of 
atomic hydrogen. The circular cross-section of the cylinder is distorted by unequal x and y scales 
into an ellipse. See text for details. 
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FIG. 4. The density distribution of Ceo molecules falling through a cylindrical condensate of 
atomic hydrogen at 45 /im below the cylindrical center. The fullerenes are released at 45 ^m above 
the trap center with initial uniform density no. The condensate's cylindrical radius is 7.5 p. 
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